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ABSTRACT 



1 Preprint online version: August 17, 2012 

o 

Context. The Sloan Digital Sky Survey (SDSS) allows us to classify galaxies using optical low-ionization emission-line diagnostic 
diagrams. A cross-correlation of the SDSS data release 7 (DR7), containing spectroscopic data, with the Very Large Array (VLA) 
survey Faint Images of the Radio Sky at Twenty-Centimeters (FIRST), makes it possible to conduct a joined multiwavelength statis- 
■ tical study of radio-optical galaxy properties on a very large number of sources. 

Aims. Our goal is to improve the study of the combined radio-optical data by investigating if there is a correlation between the radio 
luminosity at 20 cm over the luminosity of the optical Ha line (L 2 o cm I L Ha ) and line excitation ratios, where the latter provide the 
spectroscopic classification in Seyferts, low-ionization nuclear emission-line regions (LINERs) and star-forming galaxies. We search 
for a trend with z in the classification provided by classical and more recent optical emission-line diagnostic diagrams. 
Methods. We cross-match the optical sources with the FIRST radio survey in order to obtain spectroscopic information of a selected 
sample of radio-emitters with optical observed counterpart. We search for a L 2 o C mlL Ha threshold value above which the radio emit- 
ters start being classified as Active Galactic Nuclei (AGNs) rather than star-forming galaxies (SFGs). We investigate the origin of 
emission-lines by using both photoionization and shock models. 

Results. The percentage of detected AGNs (Seyferts and LINERs) or composites is much higher in the optical-radio sample than 
in the optical sample alone. We found a progressive shift of the sources towards the AGN region of the diagram with increasing 
L20 cm/Luc, with an indication for a different behavior for LINERs and Seyferts. The classification appears to slightly depend on the 
redshift. The diagnostic diagrams display a density peak in the star-forming or composite region for log(L 2 o cmlL Ha ) < 0.716, while 
the distribution in the LINER region peaks above this treshold. A comparison with photoionization and shock models shows that the 
\ large fraction of LINERs identified in our study have emission lines that may be explained by shocks. 

Conclusions. Our results indicate that it is worthwhile to further explore the radio domain, probing the physical nature of LINERs 
thanks to a combination of optical and radio information. The [Nii]/Hq' vs. equivalent width of the Ho- line (WHAN) diagram confirms 
the LINER classification for most of those that have been identified with the traditional diagnostic diagrams. The correlation between 
^ ' L20 cmlLna and optical emission line ratios suggests the nuclear origin of the emission from the most powerful radio-galaxies. 

m 
en; 

qq ■ 1. Introduction can be detected at higher redshift. Radio galaxies can also be 
C — ) , , , , , , . , , . classified according to whether they have strong high-excitation 
^ ■ Our knowledge of global galaxy properties has been recendy im- narrow . line emi ssion. The majority of FR I radio galaxies show 
^ . proved thanks »to ^lar ge-area surveys such as the Sloan Digital wgak Qr letel absent emi ssion lines. Tho se are re- 
j. . Sky Survey (SDSS; | Y rk et al lB, which have allowed for fefred tQ ag low . excitation tems (Ha rdcastle et al. 200$ and 
>■ statistical studies on large samples of galaxies. SDSS data can ^ ^ mosd found in em ^ laxigs wkh Utde 
^ . be used as a starting point to conduct multiwavelength studies . ^ formation (Ledlow & Qwen 1995; Govoni et al 2000) 
1 ° n 8 alax y Properties outside the optical regime such as the radio Tf ^ ctro ic information is available, these galax- 
Oj ; domain, where Active Galactic Nuclei (AGNs) can be detected. ieg ^ „ classified as low-ionization nuclear emission- 
Radio ga axies are active galaxies that are very luminous at ra- ^ ionf , (LINERs , Heckman M^)- Conversely, the most 
dio wavelengths. They are considered as radio-loud if they have rful high . redsh if t F R II radio galaxies have, in most of 
ratios of radio (at 5 GHz ) to optical (B- band) flux greater feajglO the caseSi strong emission i ines and are c i assifi ed as AGNs 
d Kellermann et al. | 1 989|). According to lFanaroff & Rilev UM), ^Koziel-Wierzbowska & Stasihska 201 lb w ith peculiar optical 
radio-loud galaxies can be divided into two major classes. In hol ieSi tails> brid and shells (Smith & Heckman 

Fanaroft-Riley Class I (FR I) sources, the radio emission peaks , non , ,^1 1 a. „ • * «■ 1 T< . u 1 

, J , v , , . .' , . „ . .. 1989) and bluer colors with respect to giant ellipticals. It has also 

near the galaxy nucleus and the emission from the jets fades with been estabHshed that rful FR IIs show a sfr c<jrrelation 

distance from the centre, while class II (FR II) sources present between ^ radio luminosit and their tical emission . line 
bright radio lobes. Class I sources dominate the population of luminosit fe aum et al] |i995). suggesting that both the optical 
radio-emitters at low radio power and low redshifts while more and ^ radio emission m ori inatin in the same h ical 
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powerful radio galaxies are almost exclusively FR II systems and 
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Send offprint requests to: Mariangela Vitale, e-mail: The possibility of a more general correlation between the 
vitale@phl.uni-koeln.de emission-line and radio luminosities of radio AGN has been 



1 



M. Vitale et al.: Classifying Radio Emitters from the Sloan Digital Sky Survey 



alread y explored in the past (McCart hy! 1 19931: IZirbel &~B aum 
119951) . Radio-loud AGNs are most likely to display emission 
lines in galaxies with low velocity d ispersions and at rad io lu- 
minosities greater than 10 25 WHz' 1 dKauffmann et al.ll2008h . A 
similar correlation has been observed between the emission-line 
luminosity and the i onization state of the gas for a sample of 
low-z radio galaxies (Saunders et all 19891) . where the larger val- 
ues of the emission-line luminosities have been measured for 
the more powerful radio sources, as an indication of the pres- 
ence of a strong ionizing AGN-like field. However, a problem 
arises since the most powerful radio galaxies are generally de- 
tected at higher redshifts than the less powerful radio galax- 
ies, due to selection effects. This makes it difficult to establish 
whether the correlation is between the emission-line luminosity 
and radio luminosi ty, rather than b etween emission-line lumi- 
nosity and redshift dMcCarthvll 19931) . The well-known effect of 
having an increasing number of detected powerful radio AGNs 
with increasing redshift is not only due to selection effects, but is 
is als o supported by the downsizing scenario of galaxy evolution 
(e.g. iThomas et alj|2005l: ICimatti et al.ll2006l) . Galaxies placed 
at higher redshift are more massive and host massive black holes 
that accrete producing powerful jets. These are easily detected 
in the radio domain, while low-redshift radio galaxies host low- 
mass black holes, from which only weak jets originate. The 
quasars detected in the radio but not in the optical regime might 
be heavily obscured objects. They indeed experience dust obscu- 
ration when the line of sight passes through the optically thick 
torus, which surrounds the black hole (unified model, Antonucci 
1993). 

One of the current scenarios of galaxy evolution includes the 
possibility of a smooth transition from late-type galaxies (LTGs), 
which show blue colors and ongoing star-formation, to red and 
passively evolving early-type galaxies (ETGs). This transition 
could be dri ven by the so-called AGN-feedback (iCattaneo et al.l 
120071 120091) . which is thought to be responsible for the quench- 
ing of the star-formation. After the quenching, a phase of passive 
evolution starts and the galaxies color turns into red. Some stud- 
ies show that there is indeed a correlation betwe en the age of 
galax y stellar populations and the AGN activity dTortora et al.l 
2009, Vitale et al. 2012, in prep.), with the oldest stars inhabit- 
ing the AGN-like galaxies. Moreover, AGNs are found to reside 
almost exclusively in massive galaxies with structural proper - 
ties similar to normal early-type systems (Heck man et alJ l2004). 
which are dominated by old stellar populations. 
The question of whether the AGN-feedback is also capable of 
triggering star formation in the vicinity of the black hole has still 
to be answered. The standard acc retion mode of AGNs, w hich 
is associated with quasar activity ( Sha kura & Sunvaevll 19731). is 
relate d to star-formation in the host galaxies dKauffmann et al.l 
120031 This scenario seems to be in contrast with the quenching 
of the starburst activity due to AGN-feedback. Star-formation 
and accretion around the black hole, both fueled by the same 
material, might occur with delays with respect to each other 
dWild et al.l l2010t iTadhunter et al.l 1201 ll) and come to an end 
once the gas is exhausted. A fraction of these AGNs that ac- 
crete in "quasar" mode show powerful radio jets and are there- 
fore radio-loud. 

As an observational proof to this scenario, llvezic et al .1 d2002l) 
have shown that a sample of optically unresolved radio sources 
from the FIRST survey has bluer colors compared to other SDSS 
objects and Richards et al. (2002) found that the SDSS quasar 
candidates - which are likely to have a radio counterpart - dis- 
play blue colors, especially at z > 1, indicative of ongoing star 
formation. 



In order to understand the importance of the AGN-feedback in 
the current scenarios of galaxy evolution, a study of the AGN 
populations, the properties of the host galaxies and the physical 
mechanisms that trigger the production of the emission lines or 
the radio-activity, is necessary. 

AGNs can be selected from a spectroscopic survey using 
some optical emission line ratios. Emission-line diagnostic di- 
agrams have been extensively used during the past decades 
to point out the connection between the ga laxy nuclear activ- 
ity, its m orphological type (|Ho et alJ fl997) and its evolution- 
ary stage (tHopkins et alj|2006|). Th e Baldwin-Phillips-Terlevich 
(BPT) diagnostic diagram ( Baldwin e t alj Il98ll) and its sub- 
sequent versions (|Veil 
1 20031: lLamareille et alj 



eux & Osterbrock 1987k iKewlev eta! 
2004; Dopita et alj 120021: iGroves et al 



2W4M make use of emission line ratios whose strength is a 



function of the hardness of the ionizing field of the galaxy, the 
ionization parameter U and the metallicity (see Sec. 3.4). Higher 
ratios are thought to be mainly the product of the ionization that 
arises due to accretion around the black hole (which implies the 
AGN at the center of the galaxy is in its active phase) rather than 
photoionization by hot massive OB stars. This diagnostic tech- 
nique, largely used in the optical wavelength regime, allows to 
discriminate between galaxies which show activity in their nu- 
clei and starbursts. 

Some attempts to link optical and radio properties of a large 
sample of galaxies by using combined spectroscopic and pho- 
tometric information have been already made in the past. Obric 
et al. (2006) found strong correlations between the fraction of 
detected AGNs at other wavelengths and opt ical properti e s such 
as flux, colors and emission-line strengths, llvezic etall ([2002) 
discussed the optic al and radio properties of ~ 30 000 FIRST 
dBecker et al.ll 995) sources positionally associated with a SDSS 
source by analyzing their colors, while iBest et ail (120051) com- 
pared the optical survey to both FIRST and NVSS radio surveys 
in order to derive the local radio luminosity functions of radio- 
loud AG Ns and star-formin g galaxies. In the SDSS Early Data 
Release dlvezic et al.l 120021) . ~ 70% of FIRST sources do not 
have an optical counterpart within 3". This is probably due to 
the fact that the majority of unmatched FIRST sources, detected 
down to 1 mJy sensitivity, are optically too faint to be detected 
in the SDSS images. Moreover, the fraction of quasars in the 
FIRST catalog seems to be a strong function of the radio flux, 
monotonically decreasing from bright radio sources towards the 
FIRST radio sensitivity limit. 

SDSS spectra have been used to compute the line strenghts of 
several emission lines to classify gala xies into sta rbursts and 
AGNs by using diagnostic diagrams. In llvezic et al.l d2002l) . the 
fraction of radio galaxies which are classified as AGNs rather 
than starbursts (~ 30%) is 6 times larger than the corresponding 
fraction for all the SDSS galaxies. Furthermore, the radio emis- 
sion from AGNs results to be more concentrated with respect to 
radio emission from starburst galaxies, as we expect from the 
nuclear origin of AGN emission. Radio emission is point-like 
in compact quasars detected at high redshift, while local galax- 
ies tend to have larger radio sizes. This suggests that a signifi- 
cant fraction of the radio emission either originates outside the 
nuclear region or that the radio lobes are resolved. Depending 
on the beam size, it can happen that not all the light from the 
galaxy is taken into account. The relative number of AGNs de- 
creases with the radio flux and this is consistent with the differ- 
ences in the radio luminosity functions of starburst and AGNs 
dMachalski & GodlowskfeoOOtlSadler et al.ll2002l) . 
In a multiwavelength approach, we made use of both optical - 
for the part concerning emission lines and diagnostic diagrams - 
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and radio data to conduct a statistical study on the prospects of 
identifying radio galaxies in some well-defined regions of the 
low-ionization emission lines diagnostic diagrams by using a 
combination of radio and optical properties. A comparison of 
the spectroscopic measurements with some photoionization and 
shock models is then presented to shed light on the origin of the 
emission-lines in AGNs and star-forming galaxies. 
The paper is organized as follows: in Sect. 2 we present the opti- 
cal and radio samples and their cross-matching. We describe the 
way we obtain the final sample and its completeness. In Sect. 3 
we present our results. We exploited the diagnostic diagrams for 
the cross-matched sample, looking for a trend with L20 cmll^Ha 
and redshift. We compare the spectroscopic measurements with 
photoionization and shock models. Discussion on the results is 
provided in Sect. 4. Main findings and conclusions are in Sect. 
5. 

2. Optical and Radio Samples 

2.1. Sloan Digital Sky Survey 

The Sloan Digital Sky Survey (SDSS) is a photometric and spec- 
troscopic survey that covers one-quarter of th e celestial sphere 
in th e north Galactic cap dYork et al.l l2000t [Stoughton et al. 
I2002I) . It uses a dedicated 2.5 m wide-angle optical telescope at 
Apache Point Observatory in New Mexico, United States. The 
spectra have an instrumental velocity resolution of cr ~ 65 km/s 
in the wavelength 3800 - 9200 A range. The identified galaxies 
have a median redshift of z ~ 0.1. Spectra are taken with 3" 
diameter fibers (5.7 kpc at z ~ 0.1), which makes the sample 
sensitive to aperture effects (low redshift objects are likely to 
be dominated by nuclear emission, e.g. Kewley et al. 2003). 
Several physical parameters of the galaxies in the SDSS have 
been derived and listed in publicly available catalogs. 
The Max-Planck-Institute for Astrophysics (MPA)-Johns 
Hopkins University (JHU) Data Release 7 (DR7) of spectrum 
measurements (http://www.mpa-garching.mpg.de/SDSS/DR7/) 
contains the derived galaxy properties from the MPA-JHU 
emiss ion line analysis for the SDSS DR7 (Abazaiian et al.l 
2009). It represents a significant extension in size (~ 2 times, 
from 567 486 to 922 313 objects) and a general improvement of 
the previous DR4 data set. This is both due to improvements in 
the reduction pipeline of the SDSS and in the analysis pipeline. 
In addition, new stellar population synthesis spectra for the 
stellar continuum subtraction (updated Bruz ual & Charlotl2003l) 
are used. One caveat of the MPA-JHU DR7 is that it does not 
make a differentiation between narrow and broad emission 
lines. Since the narrow component only is used in the diagnostic 
diagrams, this fact introduces a bias. An under estimation of the 
flux could lead to galaxies misclassification. lOh et al.l (120111) 
estimated that 1.3% of the galaxies in the SDSS DR7 present 
broad components in their spectra. They also showed that 
when a double Gaussian fit (for narrow and broad components 
of the same line) is performed it is possible to recover the 
missclassified objects and increase the number of sources in the 
AGN region of the diagnostic diagrams. 

2.2. Faint Images of the Radio Sky at Twenty-Centimeters 
Survey 

The Faint Ima ges of the Radio Sky at Twenty-Centimeters 
Survey (FIRST; Becker et~allll995h makes use of the Very Large 
Array (VLA) in the B-array configuration to produce a map of 
the 20 cm (1.4 GHz) sky emission with a beam size of 574 and 



a rms sensitivity of about 0.15 mJy/beam. The survey covers an 
area of about 10 000 deg 2 in the north Galactic cap, correspond- 
ing to the sky regions investigated by SDSS, and observed ~ 10 6 
sources. At the 1 mJy source detection threshold, about one 
third of the FIRST sources s how resolved structures on scales 
of 2" - 30" (llvezic et al.ll2002l) . The FIRST catalog gives infor- 
mation about the continuum flux density peak {F pea ^) and the 
integrated flux density (F mt ) at 20 cm, which allows to separate 
resolved from unresolved sources. 



2.3. Sample selection and crosscorrelation 

The SDSS DR7 spectroscopic sample makes it possible to con- 
duct a large statistical study on a large sample of galaxies (up 
to ~ 10 6 objects). With the photometric and spectroscopic infor- 
mation publicly available, the use of such diagnostic tools, as the 
emission-line diagnostic diagrams, opens the door to new scien- 
tific interpretations. A cross-correlation with the FIRST survey 
provides us with a large optical-radio sample of galaxies. 
Galaxies at very low redshift have an angular size larger than the 
size of the fiber used for observations (for the SDSS, 3"), so part 
of the emission coming from the outer regions of the galaxies 
could be misse d. In order to take in to account the aperture effect 
and following lKewlev et al.l J2003), we consider only the objects 
with z > 0.04. 

Seyfert galaxies usually have bright emission lines and they 
are all expected to be classified in the diagnostic diagrams. On 
the other hand, LINERs present much weaker emission lines. 
Therefore, by applying any cut in EW when selecting galaxies 
for our sample, we are likely to loose a part of these optically- 
weak emitters. In order not to be biased against LINERs and 
to represent as many of them as possible in the diagnostic dia- 
grams, we opted not to apply any EW cut. Instead, we included 
in our sample only galaxies whose error on the measurement of 
the EW of the lines which appear in the diagnostic diagrams is 
less than 30%. 

Our choice was supported by several tests on the cross-matched 
sample. The weakest emission line has been found to be the 
[Oi]/16300 line. If we keep our error cut and additionally apply 
a more severe cut in the measured EW of this line (EW> 3 A 
) we are left with a quite a small number of galaxies (1 808). 
The second most affected line is [Om]/I5007 (6 138 galaxies 
left), followed by Hfi (7 573 galaxies left). The possibility to 
set this condition for all the three weakest emission lines has not 
been taken into account, since it would imply working with a 
rather small sample of radio-emitters (1 725 galaxies) and poor 
statistics. The same cut on the strongest lines - Ha, [Nn]/16583 
and [Sii]/L16717, 6731 - does not exclude a significant number 
of galaxies. Moreover, in a cross-matched sample where we ap- 
ply the redshift cut only (34 733 galaxies), the mean logarith- 
mic values of the emission line ratios are higher (log [Nul/Ha = 
-0.13, log [Sii]/Ha = -0.29, log [Oi]/Ha = -0.91) compared 
to the case in which we also apply the 30% error cut on all 
the emission lines (9.594 galaxies, log [Nn]/Ha = -0.24, log 
[Sul/Ho' = -0.42, log [Oi]/Ha = -1.22), indicating in the sec- 
ond case a shift in the classification from LINERs to starburst 
galaxies. H/3 and [Om]/I5007 appear in all the diagnostic dia- 
grams. Their ratio is strongly affected by the weakness of the 
components, especially the H/3 line. This line is often superim- 
posed on an absorption component, making its strength depen- 
dent on the quality of the stellar continuum subtraction. In the 
case of a sample in which a cut in redshift and EW error has been 
applied to all the emission lines, as well the requisite EW> 3 A 
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Emission line 


bw„, 




CDD 


net 


AC) ^(S 

4UJU 




01 

U.U 1 


up 


8.12 


0.32 


0.07 


[Om]/I5007 


15.53 


0.36 


0.08 


[Oi]/I6300 


2.21 


0.22 


0.14 


[Nn]^6583 


19.42 


0.31 


0.02 


[Sn],U6717,6731 


13.30 


0.52 


0.06 



Table 1. Statistics on the EW of the lines. In the first column we indicate the name of the measured emission lines. In the second 
column the mean values of the EW of the lines in A are reported. In the fifth and sixth columns we indicate the absolute and relative 
mean errors on the EW measurements. 



for the H/3 and [Om]/I5007 lines (5 298 galaxies), the line ratios 
slightly decrease further, moving the bulk of the population to 
the SF region. This is probably due to the lack of LINERs. The 
log [Oiii]/H/? ratios increases from -0.01 of the precedent case 
to 0.12 instead, becoming more characteristic of Seyfert galax- 
ies. This leads to the conclusion that we miss a fraction of the 
weakest optical-emitters (e.g. LINERs) when we apply an EW 
error cut on all the emission li nes, but we are still le ft with a sig- 
nificant number of sources. In iKewlev et al.l d2006l) . the authors 
prefer to consider a S/N> 3 to be sure that the quality of the 
galaxy spectra is high enoug h to make accurate line measure- 
ments. Stasihs ka et alj ((2006) mentioned that this cut does not 
change the visual shape of the galaxies distribution in the diag- 
nostic diagrams, but it reduces the proportion of the objects in 
the right wing of the seagull. 

After applying redshift and EW error cut (z > 0.04, error on the 
EW measurements of all the lines which are used in the diagrams 
is less than 30%), the resulting sample of SDSS galaxies with 
emission line measurements from DR7 contains 79 919 galax- 
ies, which is about t he amount of objects that have been used by 
IKewlev et al.l d2006l) . In that case, the authors made use of a pre- 
vious SDSS data release (DR4) containing less sources, thus our 
selection criteria are more restrictive. Tab. Q] contains the statis- 
tics on the EW of the lines for the objects that have been included 
in our final cross-matched sample. 

For generating the cross-matched FIRST/SDSS sample, we used 
the matching results p rovided by the SDSS DR7 via Casjobs 
(lOMullane et al.ll20Q5h . based on a matching radius of 1". The 
resulting sample contains 37 488 radio-emitters and represents 
nearly 40% of the ~ 10 5 FIRST sources and around 4% of the 
~ 10 6 SDSS sources. SDSS spectra are available for all ob- 
jects of the matched sample. Some objects have been lost during 
the matching procedure due to the centering problems of radio 
emission with respect to the optic al emission, esp ecially in case 
of extended sources. Moreover. iBest et ail {2005) found in their 
study that radio galaxies extracted from the main spectroscopic 
sample of the SDSS reside in very massive early-type galaxies, 
with weak or undetectable optical emission lines. This further 
proves that the cross-matching process implies the loss of many 
radio emitters that do not have optical counterparts with emis- 
sion lines. 

Another considerable fraction of radio-emitters which have an 
optical counterpart is lost after applying the redshift cut (z > 
0.04) and the error cut on the EW of the lines involved in the di- 
agnostic diagrams. Our final cross-matched sample consists of 
9 594 objects, which corresponds to 25.6% of the full cross- 
matched optical-radio sample, ~ 9.6% of the original radio sam- 
ple (FIRST) and ~ 1% of all the galaxies in the MPA-JHU data 
release. 

The radio luminosity has been calculated from the integrated 
flux provided by the FIRST catalog, assuming a cosmology with 
H = 70 km &- 1 Mpc" 1 , Q M = 0.3, Q. A = 0.7. 



2.4. Completeness of the sample 

Our sample completeness is strongly related to the complete- 
ness of the parent samples (SDSS and FIRST). The SDSS 
spectroscopic sample can be considered as complete, in the 
sense that the biggest incompleteness comes from galaxies mis- 
classifications due t o mechanical spectrograph constraints (6%, 
iBlanton et alj |2003). which causes slight underrepresentation of 
high-density regions. For few sources (< 1%) the redshift can- 
not be determined or it has been derived in an incorrect way. 
Moreover, few targets (~ 1% ) are contaminated by galactic stars. 
According to lBlantorJ (12006). the mechanical constraints are re- 
lated to the impossibility of placing the fibers sufficiently close 
to each other. When two galaxies are found to have a small sepa- 
ration, only one source is chosen (independently from its magni- 
tude or surface brightness) and this does not result in a luminos- 
ity bias. Reliability and accuracy of the catalog of r adio sources 
extrac ted from the FIRST images is discussed in IWhite etaD 
d!997h . 

3. Results 

A statistical study of the spectral properties of radio emitters 
from the SDSS has been carried out starting from the determi- 
nation of the fraction of star-forming galaxies and AGNs com- 
posing our sample. For this purpose, we exploited the emission- 
line diagnostic diagrams to separate between sources with a hard 
ionized spectrum (AGNs) and those that are dominated by star- 
formation, characterized by a considerably softer ionizing field. 
Since all the emission lines which are used in the diagrams are 
placed close to each other in the spectrum, their ratios are almost 
insensitive to reddending. Therefore, we do not need to correct 
for the interstellar and galactic extinction. The luminosity of the 
Ha line has been derived after correcting the corresponding flux 
for the visual extinction by using a theoretical Ha/Hp Balmer 
ratio of 2.86. 

3.1. Diagnostic Diagrams for the SDSS and the 
cross-matched sample 

Figure [T]presents the optical emission-line diagnostic diagrams 
for the MPA-JHU and the cross-matched radio-optical samples. 
The SDSS targets are represented in red, while the radio emit- 
ters are overplotted in orange. In all the diagrams, 79 235 of 
the 79 919 (99.1%) galaxies in our MPA-JHU sub-sample are 
represented, and 9 408 of the 9 594 radio-emitters (98,1% of 
the optical-radio sample). The difference is due to the lack of 
some emission-line measurements in the parent catalog of SDSS 
galaxies and especially the [Oi] line, which is weaker than the 
other lines. 

The dashed demarcation curve in the plots of the first column of 
Fig[U([Nii] or BPT diagram) has been derived by IKewlev et al.l 
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Fig. 1. Optical emission-line diagnostic diagrams. The SDSS (DR7) targets are represented in red, while th e radio emitters from 
the cross-mat ched sample are plotted in orange. Demarcation curves in the left panels ([Nn] diagram) are bv lKewlev et alj d2001l) 
(dashed) and Kauffman n et al.l (120031) (dot-dashed); in the middle ([Sn] diagram) and right ([Oi] diagram) panels, the demarcation 
curves are bv lKewlev et a l. (2006). The top panels show the SDSS targets overplotted with the radio-emitters. The middle and botton 
panels show the density function of our data distribution. Density levels represent 700 galaxies per contour in case of SDSS targets, 
while it is 70 for the radio-emitters. 



Diagram 


Spectral Type 


SDSS DR7 


Cross-match 


[Nil] 


AGNs 
Composites 
SFGs 


11 756±108 (14.8%) 
13 856+118 (17.5%) 
53 623+232 (67.7%) 


2 968+54 (31.5%) 

3 021+55 (32.1%) 
3 419+58 (36.3%) 


[Sn] 


Seyferts 
LINERs 
SFGs 


7 646±87 (9.6%) 
4 107±64 (5.2%) 
67 482±260 (85.2%) 


1 747+42 (18.6%) 
1 012+32 (10.8%) 
6 649+81 (70.7%) 


[Oi] 


Seyferts 
LINERs 
SFGs 


10 269±101 (13.0%) 

6 543+81 (8.3%) 
62 423+250 (78.8%) 


2 180+47 (23.2%) 
1 500+39 (15.9%) 
5 728+76 (60.9%) 



Table 2. Statistics on the SDSS DR7 sample and our optical-radio sample. The table shows the number of different objects placed 
in the AGN, composite or SF region (left plots, FigJTJ or Seyfert, LINER or SF region (middle and right plots) of the diagnostic 
diagrams. The relative fraction of different spectral types is reported in parenthesis. The number of active galaxies is larger for the 
cross-matched sample. 



5 



M. Vitale et al.: Classifying Radio Emitters from the Sloan Digital Sky Survey 



(200lj) by constructing a detailed continuous starburst model 
with large realistic metallicity and ionization parameter ranges, 
in order to find an upper limit for the position of the star-forming 
galaxies in the diagram. This upper limit can be used to sep- 
arate AGNs from st ar-forming galaxies. Th e dot-dashed curve 
has been derived by Kauffm ann et al.l (120031) . who proposed an 
empirical and more conservative cut to identify starbursts by us- 
ing the large sample of emission line galaxies in the SDSS. This 
method selects many less star-forming galaxies than AGNs - ac- 
cording to the main ionizing mechanism that produces the emis- 
sion lines - compared to the Kewley's criteria. The enclosed re- 
gion between the two curves is considered to be populated by 
mixed or transitional objects. In the middle ([Sn] diagram) and 
right ([Oi] diagram ) panels, the demarcation curves are from 
iKewlev et al.l d2006l) . A further line allows us to separate Seyfert 
galaxies from LINERs. 

In Tab.|2]we present the statistics on the samples. For each diag- 
nostic diagram, we report the number of AGNs (with the distinc- 
tion between Seyferts and LINERs, when available), composites 
and star-forming galaxies for both samples. The absolute num- 
ber of objects per region is followed by the purely statistical error 
(Poisson error estimated as y/N) and the percentage with respect 
to the total number of galaxies with measured emission-lines. A 
higher relative number of radio-emitters are placed in the transi- 
tional or AGN regions compared to the full optical sample and 
the percentage of radio-emitters in the left arm of the 'seagull' 
shape defined by the SDSS targets ([Nn] diagram) is lower. The 
SDSS is strongly dominated by star-formation rather than AGN- 
like emission. Fig. Q] shows that the AGNs and the radio galax- 
ies with AGNs are drawn from a population which has higher 
metallicity than the overall SDSS sample. The lower abundance 
objects are predominantly-star forming, and populate the upper 
left-hand portions of these diagrams. The [Oi] diagram is the 
one that shows the largest number of objects classified as AGNs 
(39.1% of the cross-matched sample). 

We should take into account that the redshift we are considering 
in our study is limited due to the choice of the diagnostic dia- 
grams. Some of the emission lines are shifted out of the SDSS 
spectral range for objects with z *Z 0.4 (~ 2% of the optical-radio 
sample). 



3.2. Diagnostic Diagrams for the cross-matched sample: 
Trend with L 20 cm /L Ha 

The luminosity of the Ha line, Lff a , is conside red to be a good 
optica l star-formation rate (SFR) indicator dMoustakas et al.l 
2006). The ratio between the radio luminosity and Ln a can 
be used to compare emission from radio-components with 
the emission from young stars. We divided the cross-matched 
optical-radio sample in four bins with increasing Laocm/Lna 
containing approximately 2 350 galaxies each, in order to keep 
the number of objects per bin constant and to search for a 
threshold above which the objects start to be classified as AGNs 
(Seyfert or LINERs) rather than star-forming galaxies. 
We would expect to find a correlation between the optical and 
the radio emission, according to the hypothesis that they are 
linked to each other (e.g. L20 cm - L„ relation, where L„ is the 
luminosity in the optical) because they originate from the same 
physical process. Strong emission lines, such as [Om]/I5007, 
[On]/l3727 and H/3, are indeed believed to c o me from pow- 
erful ra dio-emitters dBaum & HeckmarJ fl989t IRawlings et"a"fl 
19891: IRawlings & Saunders! I199U iMorganti et al.l 119921. 
Zirbel & Bauml 119951: iTadhunter et all 119981: iBest & Hec kman 



120121) . These lines are also used to select AGNs from a spec- 
troscopic survey, pointing to a possible correlation between 
the AGN-detection rate and the radio luminosity of the host 
galaxies. 

We found that the peak of the distribution of the radio emitters 
shifts from the SF region of the diagrams to the compos- 
ite or AGN part (on the right-hand side of the diagrams) 
for increasing log(L20 cm/Lna) (Figs. [2] [3] @). In particu- 
lar, the upper-left panels mostly show starbursts with high 
metallicity, while the bottom-left panels display a mixed 
population and the bottom-right panels show a nearly pure 
AGN population together with some metal-rich starbursts. 
The distribution shows a peak in the LINER region for 
log(L20 cm/Lna) > 0.716, where 50% of the radio-emitters 
are classified as LINERs in the [Oi]-based diagram (36.2% 
in the [Sn]-based diagram, see Tab. [3}. In the upper-right and 
bottom-left panels (log(0.702 < T-qr cmlLjia) < 0.708 and 
log(0.708 < L 2 Qcm/LHa) < 0.716) of each figure, the Seyfert 
region appears increasingly more populated (from few per 
cent up to more than 30%), while for log(L2u cm/Lna) > 0.716 
the number of Seyfert galaxies seems to remain constant. On 
the contrary, the increase in the number of LINERs becomes 
exponential in the last bin, where it goes from the few per 
cent of the first 3 bins to half of the entire population of 
radio-emitters. All diagrams exhibit this behaviour, though the 
trend is most obvious in the [Oi] diagram, which is also the 
most sensitive to shocks. In Tab. [3] and Fig. [5] we report the 
statistics for the L20 cml^Ha bins on all the diagnostic diagrams. 
The number of Seyfert galaxies and AGNs in the BPT diagram 
increases progressively by a factor 2.7 in the first 3 bins, while it 
remains almost constant in the last bin. The number of LINERs 
drastically increases (x5) from the third to the last bin and, on 
the other hand, the number of SFGs decreases in the last bin. For 
each diagnostic diagram we indicate the number of classified 
star-forming galaxies, Seyferts and LINERs (or composites and 
AGNs in the [Nn] -based diagram) per L20 cm/^f/a bin. 

3.3. Diagnostic Diagrams for the cross-matched sample: 
Trend with redshift 

We further exploited the diagnostic diagrams as a function of 
redshift by dividing the sample in four bins containing approx- 
imately 2 350 objects each. As for the LzOcmlLna bins, this 
choice comes from the intent to keep the number of objects per 
bin constant, together with the intention to further explore the 
dependence of the optical classification with redshift. 
We found that the distribution of the radio-emitters in the diag- 
nostic diagrams slightly depends on the redshift (Figs. [6] [7] [8}. 
The number of AGNs (both LINERs and Seyferts) always in- 
creases with z. The bulk of the population remains in the SF, 
composite or SF+composite region in all z bins (see statistics on 
Tab. Hand Fig.g). The number of Seyferts, LINERs and AGNs 
in the BPT diagram is nearly double in the last bin with respect 
to the first one. The number of composites remains almost con- 
stant, while the number of SFGs decreases of a variable percent- 
age (22% in the [Sn] diagram, 34% in the [Oi] diagram and 45% 
in the [Nn] diagram). 

3.4. The [Nn]/Ha vs. equivalent width of the Ha line (WHAN) 
diagram 

Emission from some of the galaxies classified as LINERs in 
the diagnostic diagrams is nowadays thought to be triggered 
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-2.0 -1.5 -1.0 -0.5 0.0 -2.0 -1.5 -1.0 -0.5 0.0 
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Fig. 2. [Oi]-based diagnostic diagrams for the optical-radio sample. From top to bottom and from left to right, log(L20 cm /Ln a ) 
increases. The number of radio-emitters per bin is constant and equal to 2 350 ± 25. The contours represent the number density of 
the radio-emitters (20 galaxies per density contour). The bulk of the population shifts from the SF region to the composite or AGN 
region with increasing L20 The transition seems to occur for log(L20 cml^Ha) > 0.716. 



Diagram 


Spectral Type 


L 20 a„/L Ha < 0.702 


0.702 <L 20cm /L H „< 0.708 


0.708 <L 20cm /L H „< 0.716 


L20 cmlLhla > 0.716 


[Oi] 


Seyferts 
LINERs 
SFGs 


126±11 (5.4%) 
20±5 (0.8%) 
2 199±47 (93.8%) 


347+19 (14.5%) 

64+8 (2.7%) 
1 973+44 (82.8%) 


859+29 (36.5%) 
254+16 (10.8%) 
1 241+35 (52.7%) 


848+29 (36.5%) 
1 162+34 (50.0%) 
315+18(13.5%) 


[Sn] 


Seyferts 
LINERs 
SFGs 


88+9 (3.7%) 
10±3 (0.4%) 
2 247±47 (95.8%) 


261+16(10.9%) 

24+5 (1.0%) 
2 099+46 (88.0%) 


676+26 (28.7%) 
136+12(5.8%) 
1 542+39 (65.5%) 


722+27 (31.0%) 
842+29 (36.2%) 
761+28 (32.7%) 


[Nn] 


AGNs 
Composites 
SFGs 


134+12 (5.7%) 
659+26 (28.1%) 
1 552+39 (66.2%) 


346+19 (14.5%) 
820+29 (34.4%) 
1 218+35 (51.1%) 


918+30 (39.0%) 
925+30 (39.3%) 
511+23 (21.7%) 


1 570+40 (67.5%) 
617+25 (26.5%) 
138+12(5.9%) 



Table 3. Statistics on the L20 cm/Lna (in logarithmic value) bins of the optical-radio sample. The table shows the number of different 
objects placed in the Seyfert, LINER or SF region (Fig |2l31 l. or AGN, composite or SF region (FiglUi of the diagnostic diagrams. 
The relative fraction of spectral types is reported in parenthesis. 



Diagram 


Spectral Type 


0.04< z <0.07 


0.07< z <0.09 


0.09<z<0.13 


z>0.13 


[Oi] 


Seyferts 
LINERs 
SFGs 


385+20 (16.4%) 
288+17 (12.3%) 
1 675+41 (71.3%) 


461+21 (19.6%) 
314+18(13.3%) 
1 580+40 (67.1%) 


585+24 (24.9%) 
386+20 (16.5%) 
1 374+37 (58.6%) 


749+27 (31.7%) 
512+23 (21.7%) 
1 099+33 (46.6%) 


[Sn] 


Seyferts 
LINERs 
SFGs 


308+17 (13.1%) 
213+15 (9.1%) 
1 827+43 (77.8%) 


364+19 (15.5%) 
211+14 (9.0%) 
1 780+42 (75.6%) 


465+22 (19.8%) 
261+16(11.1%) 
1 619+40 (69.0%) 


610+25 (25.8%) 
327+18 (13.8%) 
1 423+38 (60.3%) 


[Nn] 


AGNs 
Composites 
SFGs 


553+23 (23.5%) 
758+27 (32.3%) 
1 037+32 (44.2%) 


620+25 (26.3%) 
744+27 (31.6%) 
991+31 (42.1%) 


790+28 (33.7%) 
740+27 (31.6%) 
815+28 (34.7%) 


1 005+32 (42.6%) 
779+28 (33.0%) 
576+24 (24.4%) 



Table 4. Statistics on the z bins of optical-radio sample. The table shows the number of different objects placed in the Seyfert, 
LINER or SF region (Fig |7l6b or AGN, composite or SF region (FigO of the diagnostic diagrams. The relative fraction of spectral 
types is reported in parenthesis. 
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Fig. 3. As in Fig|2]but for for the [Sn]-based diagnostic diagram. 
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Fig. 4. As in Fig[2]but for for the [Nn]-based diagnostic diagram. 
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Fig. 5. Number of identified spectral types in the L20 cm/Lna bins P er eacn diagnostic diagram. Star-forming galaxies are in orange, 
Seyferts in red and LINERs in blue (composites are in red and AGNs in blue in the [Nn]-diagram). The L20 cml^Ha value is indicated 
in the x-axis. 
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Fig. 6. [Oi]-based diagnostic diagrams for the optical-radio sample. From top to bottom and from left to right, z increases. The 
number of radio-emitters per bin is constant and equal to 2 350 + 5. The contours represent the number density of the radio-emitters 
(20 galaxies per level). The bulk of the population is always placed in the SF region, though the number of AGNs slightly increases 
with z. 



by post-asymptotic giant branch (post-AGB) and white dwarf 
stars, which are abundant in early-type galaxies. The radiation 
from old ste llar populations in the se galaxies, classified as 're- 
tired' (RGs, Stasih ska et al.ll2008l) . is harder than the radiation 
produced by young stars, providing higher emission-line ratios 
than those typical of star-forming regions. Because of the tight 
correlation between metallicity, Z, and ionization parameter, U, 
[Nii]/Ha can be used as an empirical measure of the gas metal- 



licity up to [Nn]/Ha = 0.4, where higher v alues become an in- 
dication for the presence of AGN activity (Ivan Zee et"ai] fl998; 
iPettini & Pa"gell2004t IStasihska et al. 2006). On the other hand, 
the Ha equivalent width, EW(Ha'), is considered as a power- 
ful star-formation indicator, although it is also related to strong 
non-thermal excitation (e.g. the one from Seyferts). Therefore, 
the [Nii]/Ha line ratio and EW (Ha) can be combined toge ther in 
a diagnostic diagram (WHAN, Cid Fernan des et al.l l2010). 



9 



M. Vitale et al.: Classifying Radio Emitters from the Sloan Digital Sky Survey 



. 0.04<z<0.07 



1.0 



Seyferts / 




-1.0 -0.5 0.0 0.5 

og [S!l]/Ha 



1.5 


0.07<z<0.09 ' 


. Seyferts / ' 


1.0 




. / 








0.5 






0.0 






0.5 




\ LINERs 

\ 


1.0 


_ SF 


\ 




-1.0 -0.5 


0.0 0.5 



log [Sl!]/Ha 



" 0.0 



0.09<z<0.13 


/' ' 

Seyferts / 




1.5 


z>0.13 


_j^_^Seyferts / / 




■ / 

■ X. / 
• " . -J / 


I 

\ 


1.0 

0.5 




-■ W3EM\ . / - 
I 4 / '' 

/' r\ 




S) 7 


log [Olll] 


0.0 








W \ LINERs 
\ 




-0.5 




W///J ''■ INHs 


_ SF 






-1.0 


_ SF 


\ 



.0 -0.5 0.0 0.5 
log [Sll]/Ha 



-1.0 -0.5 0.0 0.5 
log [Sll]/Ha 



Fig. 7. As in Fig(6]but for the [Sn]-based diagnostic diagram. 
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Fig. 8. As in Fig|6]but for the [Nn]-based diagnostic diagram. 
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Fig. 9. Number of identified spectral types in the z bins for each diagnostic diagram. Star-forming galaxies are in orange, Seyferts 
in red and LINERs in blue (composites are in red and AGNs in blue in the [Nn] -diagram). The redshift is indicated in the x-axis 




Fig. 10. WHAN diagnostic diagrams for the cross-matched (orange) and the SDSS (red) samples. The upper-left panel shows the 
radio-emitters of the SDSS superimposed on the SDSS parent sample. The lower panels represent samples where a less severe error 
cut (only on [Nn] and Ha) has been applied. In this way, a considerable number of objects have been recovered, which appear to be 
classified as wAGNs or RGs. Density levels represent 70 galaxies per contour in the middle panels, 500 galaxies per contour in the 
upper-right panel and 2 000 galaxies per contour in the botton-right panel. 



The WHAN diagram can discriminate between retired galaxies 
and true AGNs, besides having the advantage of being able to 
classify galaxies with weak [Oin] and/or H/3 lines. With this 
diagram it has been found that a large fraction of weak-line 
emitting galaxies in the SDSS have LINER-like emission with 
3 <EW(Hq-) < 6 A in case of actual LINERs (labelled in the di- 



agram as wAGN, where "w" stands for "weak") and EW(Ha) < 
3 A in case of retired galaxies dCid Fernandes et al. 
Galaxies that have both high values of EW(Hq') and [Nii]/Hq- 
are classified as strong AGNs (sAGNs). 

In Figure [TO] we show the WHAN diagram for the cross- 
matched optical-radio sample (orange) and for the SDSS parent 
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Fig. 11. WHAN diagnostic diagrams for the cross-matched (left column) and the SDSS samples (right column, both in orange) 
overplotted with AGNs selected from the [Oi] -based diagram. The upper-left plot shows the Seyferts (in red) among the radio- 
emitters, while the upper-right plot shows the Seyferts in the parent SDSS sample. The lower panels show the LINERs distribution 
(in blue). Contours number density is referred to the underlying orange distribution of each plot and it is equal to 70 galaxies per 
contour in the left panels and to 700 galaxies per contour in the right panels. 



sample (red). For the subsample used in this work, radio-emitters 
and parent population (upper panels) exhibit a different distri- 
bution, with the former having the peak of the distribution in 
the sAGN region but close to the demarcation line that separates 
it from the star-forming galaxies. In the SF region, the number 
density drops and the contours show a quite steep gradient. The 
majority of galaxies showing LINER-like emission in the classi- 
cal diagnostic diagrams are found to be true LINERs rather than 
retired galaxies, because they are mostly placed in the wAGN 
zone. On the other hand, in the upper-right panel the distribution 
appears to be more spread across the SF region, with a lower rel- 
ative number of true AGNs (both strong and weak) and retired 
galaxies. 

The lower panels represent samples where a different error cut 
(30% on the [Nn] and Ha lines only) has been applied. In 
this way, a considerable number (12 419 from the cross-match, 
398 272 for the SDSS sample) of objects with weak [Ora], [Oi] 
and Hy6 lines have been recovered. They appear to be mostly 
classified as wAGNs or RGs (as shown by the density con- 
tours). Though we miss many passive galaxies in the lower part 
of the upper WHAN diagrams, we found that the shape of the 
upper part of the distributions is almost preserved for both the 
radio-emitters and the SDSS targets. The comparison with the 
lower panels suggestes, as already mentioned in Sec. 2.3, that 
the galaxy classification is strongly influenced by the error cut 



we apply on the samples. 

Figure fTTI shows the WHAN diagnostic diagrams for the cross- 
matched and the SDSS samples (both underlying in orange in 
the plots) overplotted with AGNs selected from the [Oi] -based 
diagram. We chose the latter as the one diagram that can classify 
the highest number of Seyferts and LINERs (see Tab. [3]). The 
Seyferts distributions peaks in the sAGN region of the WHAN 
diagram, with only few galaxies falling in the SF or wAGN 
regions. Moreover, it appears like most of the Seyferts do not 
lie close to the demarcation curve between SF and sAGN. The 
LINERs distributions look much more spread across the sAGN 
region instead, while they are supposed to be confined in the 
wAGN (in case of true LINERs) or in the RG region (in case of 
fake LINERs). Some LINERs appear in the SF region as well. 
Interestingly, the latter are found to have low [Oin]/Hy6, while the 
LINERs that fall in the RG region of the diagram have an aver- 
age higher [Oin]/Hy6. Fig.[T2]shows WHAN diagnostic diagrams 
for the cross-matched sample only overplotted with Seyferts and 
LINERs selected from the [Oi] -based diagram. Each plot rep- 
resents a different L20 cm/Lna bin (as in Figs. |2 [3] |4j. While 
the distribution of Seyferts always peaks in the sAGN region 
and only the number of sources increases, LINERs are mostly 
present in the sAGN region in the first three bins and they ap- 
pear to be mostly located in the wAGN region only in the very 
last L20 cm/Lfja bin. The radio strong (large L20 cm/Affa) sources 
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Fig. 12. WHAN diagnostic diagrams for the cross-matched sample (orange) overplotted with Seyferts (red) and LINERs (blue) 
selected from the [Oi]-based diagram. Each plot represents, from the top-left to the botton-right, a different L20 C mlLH a bin (Figs. 
l2l E] |4|. Contours number density is referred to the underlying orange distribution of each plot and it is equal to 70 galaxies per 
contour. 



deviate from the RG and passive regions, thus we can say that the 
radio-selection helps removing RGs from the sample. On the one 
hand, this leads to a cleaner pure AGN sample whose LINERs 
may well be dominated by shocks. 



3.5. Comparison with models 

With the exploitation of optical emission-line diagnostic dia- 
grams, both classical and WHAN, we stressed the presence of 
a higher concentration of LINERs with increasing L20 cm/Lfia- 
The point whether LINERs are pure AGNs - i.e. the main ion- 
izing mechanism in the host galaxy is accretion onto the super- 
massive black hole - is still matter of debate. There are several 
processes that could account for LINER emission: photoioniza- 
tion by a non-stellar source such as an AGN accreting at a low 
rate dHo et al.lll99l iGroves et alj|2004at lEracleous et alfe OlO) 
and io nization by fast large scale shocks (Dopita & Sutherland 
Il995l) produced either by AGN radio-jets, galactic super- winds 
dLehnert & Heckman|[l996t iLutz et alj|1999t ICecil et alj|2001l) 
or starburst-driven super- winds dHeckman et al.l Il990h . It has 
also been found that hot post-asymptotic giant branch (post- 
AGB) stars and white dwarfs can be responsible fo r the LINER- 
like e mission in a considerable fraction of galaxies dBinette et alJ 
G221. 

In order to further explore this topic, we decided to investi- 
gate the physical mechanisms from which the observed opti- 



cal emission-lines are originating. According to Dopita & Evans 
(1986), the emission lines in the spectrum of individual HII re- 
gions are determined by three main physical parameters, namely 
the temperature of the stars, T, and thus their age, the metallic - 
ity and the ionization parameter, U (usually defined as the ra- 
tio of the mean ionizing photon flux to the mean atom density). 
Using available photoionizati on and shocks models that can ac- 
count for the line emis sion dGroves et alJ l2004ah [Dopita et al. 
2006; lAllen etal] |2008), it is possible to unveil the dominant 
physical processes taking place in the galactic nuclei and their 
vicinity. For this purpose, we used the [Snl/Ha vs. [Nnl/Ha 
and [Oil/Ha vs. [Nul/Ha diagrams. This choice has been made 
by considering that these emission line ratio s do not depend 
strongly on the age of the stellar population dMov et alJ [200 1 ; 
IStasinska & Izotovl2 003; Sta sinska etaT]|2006l) . 
The left panels in FigQ~3]shows the models compared to the op- 
tical sample and the cross-matched optical-FIRST sample. We 
found that the position of most objects can be explained by mod- 
els with Z = 2Z . The star-formation models are represented by 
green lines. Solid lines indicate different 

logU values, ranging from (bright green) to -6 (light green- 
blue line). Dashed lines indicate various a ges, from 0.5 Myr s 
(bright green) to 3.0 Myrs (light green-blue) (Do pita et alj2 006). 
Models that account for AGN emission, parametrized by 
a smooth a nd featureless p ower law or br oken power law 
dKos]dl978l:IStasinskall984tlOsterbroc]dl989t lLa or et all 19971: 
iPrieto & Viega s 2000; S azonov et al. 2004), are mainly charac- 
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Fig. 13. Data comparison with models. Galaxies from the MPA-JHU sub-sample are in black, while the optical-radio sample is in 
red (left panels). In the middle panels, Seyferts of the optical-radio sample selected via [Oi]-based diagram are in red . In the right 
panel s, LINERs selected with the same diagram are in blue. All models have Z = 2Z e . The star-formation models (Dopit a et al.l 
2006) are r epresented by gre en lines; the dusty AGN models dGroves et al .1 l2004ah are in yellow; the models for shocks without 
precursors dAllen et al]l2008l) are in blue. In the middle panels, Seyferts of the optical-radio sample selected via [Oi] -based diagram 
are in red. In the right panels, LINERs selected with the same diagram are in blue. 



terized by the ionization parameter. Dust is related to the ion- 
ization parameter by the amount of ionizing photons that can be 
scattered. It is thought to be a fundamental component of AGN 
tori; th ere is also evidence for dust in the narrow -line region 
(NLR; iTomono et alj 120011; Ik adomsk i et all 120031) . The dusty 
AGN model from Groves et al. I (I2004ah . represented in yellow 



in Fig.[T3] considers a ionization parameter ranging from 10 to 
10~ 4 and spectral index a = -2.0. 

The m odels for shocks without precursors are from lAllen et al] 
(2008). The authors consider only fast shocks, where the ion- 
izing radiation generated by the hot gas behin d the shock front 
creates a strong photoionizing radiation field (ISutherland et alJ 
[1991 iDopita & Sutherland|[l99l fl996h . The flux of the ioniz- 
ing radiation emitted by the shock, and so the ionization param- 
eter, increases with the velocity of the shock front (oc y 3 ). The 
magnetic fields in the gas behind the shock front act to limit the 
compression through the shock itself. This effect is treated with 
the introduction of the magnetic parameter (B/n e 1 ^ 2 ) where the 
electron density (n c ) is equal to 1 cm~ 3 . The different values of 
the magnetic parameter range from B/n e ^ 2 = (solid lighter 
blue line) to 4 /uG cm 3 ^ 2 (solid darker blue line). Velocities go 
from 100 km s _I (dashed lighter blue line) to 300 km s^ 1 (dashed 
darker blue line). 

From the first row of plots in Fig.[T3l([Siil/Hg vs. [Nn]/Ha di- 
agrams) it appears that we are only partially tracing the star- 
formation in the cross-matched sample, since the left-hand side 
of the diagrams, where the green lines indicating photoionization 
by hot massive stars are located, does not represent all the radio- 
emitters. The star-forming region is indeed mostly populated by 
sources th at do not have strong detected radio emission, thus the 
models bv lDopita et al.l d2006l) represent the SDSS parent sam- 
ple better. The higher emission-line ratios of the rest of the radio- 



emitters could be explained by shocks and dusty AGNs models. 
The blue lines that extend from the central part of the diagram to 
the upper right corner represent quite well the full lenght of the 
distribution. In the upper-middle and upper-right panels of Fig. 
[T3l the Seyferts and the LINERs of the cross-matched optical- 
radio sample selected by the [Oi] -based diagram are represented. 
The emission-line ratios of the Seyferts cannot be fully explained 
by the models for shocks and dusty-AGNs that we considered 
in this work. On the other hand, LINERs are mostly placed in 
the region of the diagram which is eclosed by the sequences 
for fast shocks and dusty-AGNs, with only few sources lying 
in the region explained by photoionization by young hot stars. 
The latter are found to have a rather low [Om]/Hj6 ratio. In gen- 
eral, LINERs shift to the upper-left corner of the diagram for 
increasing [Oin]/Hj6. However, many galaxies of different spec- 
tral type are not explained by these models. The second row of 
Fig. fT3l ([Oil /Ha vs. [Nnl/Ha diagrams) shows that either the 
cross-matched nor the optical samples are well represented by 
models for photoionization by hot massive stars. The Seyferts 
and LINERs are instead much more well confined within the 
fast shocks and dusty AGNs models with respect to the [Sn]/Ha 
vs. [Nn]/Ha diagram. 



4. Discussion 

4. 1 . Classical diagnostic diagrams 

With a cross-matching of optical (SDSS) and radio (FIRST) 
data, we obtained a sample of radio-emitters with optical 
counterpart and thus optical spectral information. By using di- 
agnostic diagrams, we were able to separate between AGNs and 
starburst galaxies, studying the dependence of this classification 
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on quantities such as L20 cm/Lna or z. 



4.1 .1 . General trend of the full optical and optical-radio 
samples 

In our study of the optical-radio sample, we found that a rela- 
tively large number of radio emitters tend to be placed in the 
transitional or AGN region of the diagnostic diagrams (Fig. [T| 
Tab. 0, suggesting a strong connection between the detection 
of radio emission/features and the AGN activity. The percent- 
age of detected AGNs increases from the [Sn] (29,4%) to the 
[Nil] (31,5%) and [Oi] (39,1%) diagrams, partly due to the 
higher sensitivity of [Oi] to shocks. Our MPA-JHU subsample of 
SDSS galaxies populates a larger area of the SF region instead. 
The percentage of star-forming galaxies is maximum in the [Sn] 
(85, 2%) diagram, the fraction of pure SFs (without composites) 
is 67,7% in the [Nn] diagram and 78, 8% in the [Oi] diagram. 
This is mainly due to the selection of more starburst galaxies 
than in the radio-sample, where the star-formation can be ob- 
served in the optical but not in the radio, below the 1 mJy sensi- 
tivity threshold. 

The galaxy radio-emissio n seems to be either a function of opti- 
cal luminosity or redshift dHooper et alj!995h . The high number 
of AGNs among the population of radi o-emitters has been al- 
ready pointed out by other authors (e.g. [Buttiglione et al.f 2010) 
but never extensively investigated with a set of classical and new 
diagnostic diagrams and by using a large sample of galaxies. 

4.1 .2. Radio-emitters: trend with L 2 o cm/L Ha 

We divided our sample in four equally populated bins with dif- 
ferent ranges of log L20 cm/LHa- The significant variation in the 
number of galaxies justifies our bins choice, although we are 
aware that the errors on the luminosity could affect the distribu- 
tion of the sources within the bins. Although we expect a dif- 
ferent behaviour for FRI and FRII, with the latter being more 
abundant in the highest L20 cm/LHa bin, where the radio lumi- 
nosity is higher, it is not possible to prove this with FIRST data 
alone. Due to the large beam size of the FIRST radio survey, 
most of the galaxies are not structurally resolved, thus a distinc- 
tion between jet-dominated and lobe-dominated emitters is not 
made. 

The optical classification displays a peak in the star-forming or 
composite region of the diagnostic diagrams for relatively low 
L20 cm/Lna values (L20 cm/LHa < 0.716), while the distribution 
appears to have a peak in the LINER region when above this 
threshold. We observe in all the diagrams a progressive shift of 
the sources towards the AGN region with increasing L20 cm/Lna, 
indicating a change in the galaxy properties (Figs. [2] [4] and 
[3]). We think that the observed shift of the galaxy distribution, 
from the SF region of the diagnostic diagrams to the Seyfert 
region first and LINER region eventually, could be determined 
by two factors. First, the progressive increase of the radio lumi- 
nosity peak distribution of about one order of magnitude (from 
log L20 cm ~ 29.2 to log L20 cm ~ 30.2) has the effect of classify- 
ing many more objects as AGNs in general and Seyferts in par- 
ticular. In this respect, the larger number of AGNs suggests the 
radio nature of the nuclear and extended emission. Possibly, we 
are looking at low-luminosity AGNs producing radio-emission. 
Second, the different Ln a distribution of Seyferts and LINERs, 
with the LINERs always peaking for a particular bin at lower 
L^ values with respect to Seyferts, assumes its maximum sepa- 



ration in the last Lqa cm/Lna bin. Here Ln a of the LINERs peaks 
at its lowest value, 10 41 erg/s, and Ln a of the Seyferts at its high- 
est value, 10 42 erg/s), thus the high L20 cm/Ln a of the last bin is 
the reason for the appearance of a very large number of LINERs. 
The [Oi] diagram (Fig. [2j is the one that represents best the ob- 
served increase in the number of radio emitters in the LINER 
region (from 1% to 50%, see Tab. [3]) for increasing value of 
£20 cm/Lna- This could be due to the [Oi]/Ho' ratio, which is par- 
ticularly sensitive to shocks. However, the use of the diagnostic 
diagrams might present some problems. It is well established 
that the activity from the AGN and the star-formation inside the 
host are related, as they depend on the presence of fuelling mate- 
rial (Kauffmann et al. 2003a). Galaxies that do not present strong 
nuclear radio activity could be detected via radio emission aris- 
ing from supernovae explosions, which can contribute signifi- 
cantly to the emission at low radio-levels. In that case, it is likely 
that these AGNs show themselves as starbursts (low tale of the 
radio luminosity distribution) or LINERs (high radio emission) 
in the classification provided by the diagnostic diagrams. 

4.1.3. Radio-emitters: trend with z 

The finding of more AGNs at high redshift has been already no- 
ticed in the past and could be related to selection effects. At 
high redshift, we are likely to select powerful emitters such as 
Seyferts, while LINERs seem to be traced by radio-emission 
coming from shocks, so they are easier to find at higher radio 
luminosities. We explored the dependence of the optical galaxy 
classification on the redshift by dividing our sample of radio- 
emitters in four different redshift bins. The classification ap- 
pears to slightly depend on the redshift. The number of classified 
AGNs increases almost by a factor 2 from the first to the forth 
redshift bin. This is true for both kinds of AGN, with a small 
difference between LINERs (the average factor is equal to 1.65, 
based on the statistics of the [Oi] and the [Sn] diagrams, Tab. [3]) 
and Seyferts (1.95). This suggests that the AGN detection rate 
could be biased against the higher number of powerful sources, 
such as Seyferts, which are detected at higher redshifts. At the 
same time, it does not exclude AGNs and high redshift to be a in 
a genuine relation. 

4.2. Recent diagnostics 

The WHAN diagnostic diagram can help in sheding light on the 
nature of the emission-lines in the LINER-like spectra of galax- 
ies selected with the classical diagnostic diagrams. Fig.[TTIshows 
that Seyferts (selected by using the traditional [Oi] -based diag- 
nostic diagram) are mostly placed in the sAGN region and only 
few of them are located in the wAGN region. Their number den- 
sity drops at log [Nn]/Ho' < -0.3, before reaching the left bound- 
ary with the SF region. LINERs mostly appear in the wAGN re- 
gion, but they are spread in the RG and sAGN regions as well. 
One of the possible reasons for the misclassifications of these 
objects could be that we selected AGNs from the [Oi]-based di- 
agram. A full consistency between the latter and the WHAN di- 
agram cannot be achieved. 

Fig- [12 as well as Figs. [2] |4]and[3] display how the separation 
between Seyferts and LINERs becomes more pronounced in the 
last L20 cm/Lna bin, where many LINERs are identified, indicat- 
ing a general agreement between the classification provided by 
the different diagnostic diagrams. 

We found the use of the WHAN diagram complementary to the 
one of the classical diagnostic diagrams, as the former is the 
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only current optical diagnostic tool that can further distinguish 
between objects with LINER-like emission. In particular, the 
radio-sources of our sample are more likely to be identified as 
true LINERs (wAGN) than as passive galaxies, suggesting that 
photoionization by old stars possibly gives only a minor contri- 
bution to line-emission. 

4.3. Photoionization and shock models for the [Sn]/Ha vs. 
[Nn]/Ha and [Oi]/Ha vs. [Nn]/Ha diagrams 

The presence of a large number of LINER-classified galaxies 
in the classical diagnostic diagrams does not imply a high de- 
tection rate of central engine-dominated objects. In particular, 
for our radio-selected sources, the observed high fraction of 
LINERs can be due, for example, to the detection of radio jets 
as well as SNe-driven winds. A large fraction of radio-emitters 
are found to be AGNs (see Tab. |2), thus their emission-line ra- 
tios are expected to be at least partially explained by shock and 
dusty AGN models, and less by star-formation. The star-forming 
region is, indeed, mostly populated by sources that have a shal- 
low r adio emission. Mode ls for photoionization by hot young 
stars (Stasihska et al. 2006) can nicely explain the star-forming 
region (left wing of the seagull) in the classical diagnostic dia- 
grams. Regarding the right wing of the distributions, where the 
AGNs are located, a comparison with photoionization models 
(Stasihska et al. 2008) shows that radiation from old metal-rich 
stellar populations can explain the emission-line ratios which are 
typical of LINERs. On the other hand, these models do not work 
well for Seyfert, which could be explained by the presence of 
other mechanisms to produce the observed emission-lines. 
In Fig. [13] ([Sn]/Ha vs. |"Nnj/Ha d iagram , first row) the star- 
formation models by Dopit aetail (120061) represent the SDSS 
parent sample better than the cross-matched radio sample. 
Although the bulk of the radio-emitters is placed in the region 
of the diagram that is enclosed by the models for star-formation, 
shock and dusty-AGN models could account for most of the 
cross-matched sources that are found to be AGNs, and LINERs 
in particular (upper-middle and upper-right panels). However, 
many radio-emitters are placed in regions of the [Sn]/Ha vs. 
[Nn]/Har diagram that the models we chose cannot predict. This 
is especially the case for Seyferts. This points to the possibil- 
ity for these objects to have emission-lines that can be exlained 
by multiple mechanisms acting at the same time, e.g. to pro- 
duce the higher [Nnj/Ha ratio that no model can predict. The 
[Oi]/Hor vs. [Nnl/Ha diagram (Fig. [13] second row) does not 
show any Seyfert or LINER of the radio-selected sample in the 
region eclosed by the star-formation models, but it shows more 
clearly that the AGNs (most of the Seyferts and almost all the 
LINERs) are best represented by shock and dusty AGN models. 
These diagrams show that, as mentioned in Sect. 3.2, the radio 
sources are drawn from a metal rich population, whether star- 
burst, LINER or Seyfert-dominated. 

However, this is not sufficient to disentangle the complicated 
"nature of LINER" issue. What we can say is that some LINERs 
are triggered by mechanisms responsible for shocks, thus they 
are linked to radio-emission and easily identifiable with increas- 
ing L20 cm/Lno values. We acknowledge that these evidence are 
just a step to unveil the nature of LINERs and further studies are 
required to provide a thorough investigation. 



5. Conclusions 

Our main findings and conclusions are: 



- The requirement of detected radio emission predominantly 
selects active galaxies, thus a considerable amount of AGNs 
and metal-rich starburst galaxies where the radio-emission 
has mainly a stellar origin. Radio emission correlates with 
optical emission-line ratios, thus a cross-match allow us to 
classify and study radio-emitters according to their optical 
spectral identification. 

- Emission-line diagnostic diagrams show that many radio- 
emitters are classified as AGNs with increasing L20 cm/Lna- 
In particular, the increase of L20 cm selects powerful radio- 
emitters such as Seyferts, while the decrease of Ln a strongly 
contributes to the selection of a considerable fraction of 
galaxies with LINER-like emission. 

- The higher number of identified AGNs with increasing val- 
ues of z can indicate a true correlation between these quan- 
tities as well as a bias, related to the selection of the most 
powerful radio-emitters. 

- Emission-line diagnostic diagrams are an useful tool to clas- 
sify sources according to the main ionizing mechanism pro- 
ducing emission-lines, though the classification depends on 
the choice of the chemical species that is used in the diagram. 
In general, we found a good agreement between the classi- 
fications given by the classical diagnostic diagrams and the 
[Nii]/Ha vs. equivalent width of the Ha line (WHAN) di- 
agram. Diagnostic diagrams can only give hints on the na- 
ture of the observed emission-lines, making a complemen- 
tary comparison between data and models necessary. 

- While the star-forming sequence in the diagnostic dia- 
grams can be successfully fit by photoionization models, 
the AGN region seems to collect objects whose observed 
emission-lines are due to different processes. Most of the 
radio-emitters of our sample, which are mainly classified 
as LINERs at high L20 cmlLn a values, have emission-lines 
whose ratio can be explained by fast-shock and dusty-AGN 
models. Shocks, closely linked to the presence of radio emis- 
sion, can be produced both by stars and AGNs, thus unveil- 
ing the nature of LINERs requires a more detailed study. 

- By using diagnostic diagrams, it is possible to select popula- 
tions of LINER-like objects and further distinguish between 
the 'true-AGNs' and the galaxies whose emission is pro- 
duced by old stars. For our radio-selected sample, LINERs 
are more like true AGNs than retired galaxies. Resolving the 
spatial radio structure of these groups of LINER-like objects, 
as well as Seyferts that are particularly bright in the radio, 
would help us to shed light on the physical mechanisms that 
are responsible for the observed radio luminosity and optical 
emission-line ratios. 
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